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We consider a superintegrable Hamiltonian system in a two-dimensional 
space with a scalar potential that allows one quadratic and one cubic 
integral of motion. We construct the most general cubic algebra and we 
present specific realizations. We use them to calculate the energy spectrum. 
All classical and quantum superintegrable potentials separable in Cartesian 
coordinates with a third order integral are known. The general formalism is 
applied to quantum reducible and irreducible rational potentials separable 
in Cartesian coordinates in E2. We also discuss these potentials from the 
point of view of supersymmetric and PT-symmetric quantum mechanics. 

1 Introduction 

In classical mechanics a Hamiltonian system with Hamiltonian H and 
integrals of motion Xa 

H = ^QikPiPk + V{x,p), Xa = faix,p), a=l,...,n-l , (1.1) 

is called completely integrable (or Liouville integrable) if it allows n 
integrals of motion (including the Hamiltonian) that are well defined 
functions on phase space, are in involution {H,Xa}p = 0, {Xa,Xi,}p = 0, 
a,b=l,...,n-l and are functionally independent ({, }p is a Poisson bracket). 
A system is superintegrable if it is integrable and allows further integrals of 
motion Yb{x,p), {H,Yb}p = 0, b=n,n+l,...,n+k that are also well defined 
functions on phase space and the integrals! i/, Xi, ■••!^n-i,^) are 
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functionally independent. A system is maximally superintegrable if the set 
contains 2n-l functions, quasi-maximally superintegrable if it contains 2n-2 
and minimally superintegrable if it contains n+1 such integrals. The 
integrals Yi, are not required to be in evolution with Xi,...Xn-i, nor with 
each other. The same definitions apply in quantum mechanics but 
{H, Xa,Yi,} are well defined quantum mechanical operators, assumed to 
form an algebraically independent set. 

Superintegrable systems appear in many domains of physics such quantum 
chemistry, condensed matter and nuclear physics. The most well known 
examples of (maximally) superintegrable systems are the Kepler-Coulomb 
[1,2] system V{x) = ^ and the harmonic oscillator V{x) = ar^ [3,4]. A 
systematic search for superintegrable systems in two-dimensional Euclidean 
space E2 was started some time ago [5,6]. In 1935 J. Drach published two 
articles on two-dimensional Hamiltonian systems with third order integrals 
of motion and found 10 such integrable classical potentials in complex 
Euclidean space i?2(C) [7]. A systematic study of superintegrable classical 
and quantum system with a third order integral is more recent [8,9]. All 
classical and quantum potentials with a second and a third order integral of 
motion that separate in cartesian coordinates in the two-dimensional 
Euclidean space were found in Ref 9. There are 21 quantum potentials and 
8 classical potentials. 

The classical potentials were studied earlier [10]. In all 8 cases of 
superintegrable systems, separating in Cartesian coordinates and allowing a 
third order integral of motion, the integrals of motion generate a cubic 
Poisson algebra. In many cases this polynomial algebra is reducible, that is 
it is a consequence of the existence of a simpler algebraic structure. We 
have also studied trajectories and have shown that bounded trajectories are 
always closed for these superintegrable potentials. 
The quantum case is much richer : 21 superintegrable cases of the 
considered type exist, 13 of them irreducible. In this context we call a 
potential, or a Hamiltonian "reducible" if the third order integral is the 
commutator (or Poisson commutator) of two second order integrals. The 
potentials are expressed in terms of rational functions in 6 cases, elliptic 
functions in 2 cases and Painleve transcendents [11] Pi,Pii and Pjv in 5 
cases. 

The three reducible cases are 

V=^{x^ + y^), V=^{x' + y^) + ^ + ^, V = ^{Ax^ + y^) + ^ + cx. 
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The irreducible potentials with rational function are 



Potential l.V = h^[^ + + (^] 

Potential 3. F = ^{Qx^ + y"^ 

'A{^' + y') + ^ + -. 

Potential Q.V = h^[M^^ + V^) + + + — 



Potential 2. 1/ = ^(Qx^ + y^) 

Potential A.V = h^^^^^^ ^ ^— ^ 
Potential 5. F = ^^(^[(x^ + ^z^) + ^ + + 



(j/+a)2 (y—a)^ (x+a)^ {x—a)^ 



It is well know that in quantum mechanics the operators commuting with 
the Hamiltonian, form an o(4) algebra for the hydrogen atom [3,4] and a 
u(3) algebra for the harmonic oscillator. We can obtain from the algebra 
the energy spectrum. In many cases the algebra is no longer a Lie algebra 
and many examples of polynomial algebras were obtained in quantum 
mechanics [12,13,14,15,16,17,18,19,20,21,22]. C.Daskaloyannis studied the 
case of the quadratic Poisson algebras of two-dimensional classical 
superintegrable systems and quadratic (associative) algebras of quantum 
superintegrable systems [17]. He shows how the quadratic algebras provide 
a method to obtain the energy spectrum. He uses realizations in terms of 
deformed oscillator algebras [18]. Potentials with a third order integral can 
be investiged using these techniques. 

Supersymmetry was originally introduced in the context of grand 
unification theory in elementary particle physics in terms of quantum field 
theory involving a symmetry between bosons and fermions [23]. So far there 
is no experimental evidence of SUSY particles. At our energies we can 
distinguish bosons and fermions and this symmetry should appear as a 
broken symmetry. Supersymmetric quantum mechanics (SUSYQM) was 
introduced by E.Witten [24] as a toy model to study supersymmetry 
breaking. This method is related to earlier investigation of spectral 
properties of Sturm-Liouville diff'erential operators by G.Darboux [25] and 
T.F.Moutard [26] in the 19th century. SUSYQM is also related to the 
factorization method that was used by E.Schrodinger in the context of the 
quantum harmonic oscillator [27]. The factorization method was 
investigated more systematically later by L.Infeld and T.F.Hull [28]. 
SUSYQM is now an independent field with applications to atomic, nuclear. 
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condensed matter, statistical physics and quantum mechanics [29]. The 
relation with exactly solvable potentials has been discussed [30,31] and also 
with superintegrable potentials and quadratic algebras [32]. 
This paper is organized in the following way. In Section 2 we give the 
general form of the cubic algebra for two-dimensional systems with a 
quadratic and a cubic operator that commute with the Hamiltonian. We 
give a realization of the cubic algebra in terms of parafermionic oscillator 
algebras. We study the finite dimensional representations of the cubic 
algebra. In Section 3 we apply this method to the case of irreducible 
potentials separable in Cartesian coordinates in E2 with a third order 
integral. In Section 4 we investigate the irreducible potentials from the 
point of view of supersymmetric quantum mechanics. In Section 5 we give 
the generating spectrum algebra of the irreducible Potential 1. In Section 6 
we investigate the complexification of the irreducible Potential 1. (All other 
cases can we obtained from Potential 1). 



2 Cubic and parafermionic algebras 

We consider a quantum superintegrable system with a quadratic 
Hamiltonian and one second order and one third order integral of motion 



H = a{quq2)Pl + 2%i, ^2)^1^2 + c(gi, gs)^! + V{qi, gs) (2.1) 

A = d{qi, q2)P^+2e{qi, q2)PiP2+f{qi, q2)Pi+g{qi. q2)Pi+Hqi, q2)P2+Q{qi, q2) 
B = u{qi,q2)P! + 3v{q,,q2)P^P2 + 3w{qi,q2)PiP^ + x{qi,q2)P^ + j{quq2)P? 
+2k{qi, q2)PiP2 + l{qi, 92)^^ + m(gi, 52) A + n{q^, q2)P2 + S{qi, ga) , 



We assume that our integrals close in a cubic algebra. This is the quantum 
version of the cubic Poisson algebra obtained earlier [10] and the cubic 
generalization of the quadratic algebra studied by C.Daskaloyannis [17]. 



with 



Pi 



= —ihdi, P2 = —ihd2 
[H,A] = [H,B] = Q 



(2.2) 
(2.3) 
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The most general form of such an algebra is 



[A,B] = C 
[A, C] = aA^ + f3{A, B} + -/A + 6B + e 
[B, C] = ijA^ + uA^ + pB^ + a{A, B} + ^A + r]B + C 



(2.4) 



where {, } denotes an anticommutator. The coefficients a , /? and /i are 
constants, but the other ones can be polynomials in the Hamiltonian H. 
The degrees of these polynomials are dictated by the fact that H and A are 
second order polynomials in the momenta and B is a third order one. Hence 
C can be a fourth order polynomial. The Jacobi identity 
[A, [B, C]] = [B, [A, C]] implies p = —(3 , a = —a and r] = —7. We obtain 



[B,C]= nA^ + uA^ - /3B^ -a{A,B} + CA--fB + C . (2.5c) 



For the polynomials on the left and right sides of Eq.(2.4) and Eq.(2.5) to 
have the same degree we must have 



where ao, ... , Cs ^ire constants. 

The Casimir operator of a polynomial algebra is an operator that commutes 
with all elements of the algebra : 



[A,B] = C 
[A, C] = aA^ + I3{A, B} + -iA + 5B + e 



(2.5a) 
(2.5b) 



a = ao, P = Po, /i = ^J'0 
'l = 7o + liH, 6 = So + 6iH, e = £0 + eiH + e2H'^ 
iy = iyo + i^iH, ^ = eo + eii^ + 6^' 
C = Co + Ci^ + C2^' + C3i^', 



(2.6) 



[K,A] = [K,B] = [K,C] = 



(2.7) 



and this implies 



K = C^- a{A\ B} - p{A, B^} + (a/3 - j){A, B} + (/S^ - 6)B^ 
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(+/57 -2e)B + + ^ (z. + + (-1/./?^ + ^ + ^ + + ^) (2.8) 

+ (-^/u/35+^ + a7 + 2C)^ • 

Ultimatly, the Casimir operator will be a function of the Hamiltonian alone. 
We construct a realization of the cubic algebra in terms of a deformed 
oscillator algebra [17,18] A^} which satisfies the relation 

[A^, 6*] = h\ [iV, h\ = -b, b% = ^{N), bb' = ^{N + 1) . (2.9) 

$(A^) is called the "structure function". Following C.Daskaloyannis [17] we 
request $(A^) to be a real function and impose $(0) = and $(A^) > for 
> 0. We construct a Fock type representation for the deformed oscillator 
algebra with a Fock basis \n > , n=0,l,2... satisfying 

N\n >= n\n >, b*|n >= ^/^NTT)\n + 1 >, (2.10) 
6|0>=0, b\n >= ^/^N)\n - 1 > . (2.11) 

To obtain a finite-dimensional representation we request $(p + 1) = 0. 
Let us show that there exists a realization of the form : 

A = A{N), B = b{N) + b'p{N) + p{N)b . (2.12) 

The functions A(N) , b(N) et p{N) will be determined by the cubic algebra. 
We have by Eq.(2.5.a) 

C = [A, B]=b' A A{N)p{N) - p{N) A A{N)b, (2.13) 

where AA{N) is define to be AA{N) = A{N + 1) - A{N). When we insert 
Eq.(2.12) into Eq.(2.5b) we obtain two equations that allow us to determine 
A(N) and b(N) 

A A{Ny = -f{A{N + 1) + A{N)) + e (2.14) 
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aA{Nf + 2(3A{N)h{N) + -iA{N) + 5h{N) + e = 



Two distinct possibilities occur. 

Case 1 : /? 7^ 0. We find the following solution 



A{N) = ^{{N + uf--^-p (2.15) 



KN) = ^iiN + uf _ 1) + «^ - - 27^/^ + 4/?'^ 



The constant u will be determined below using the fact that we require that 
the deformed oscillator algebras should be nilpotent. Eq.(2.5c) gives us 

2$(iV + l)(AA(iV) + |)p(iV) - 2$(A^)(AA(A^ - 1) - |)p(iV - 1) (2.16) 
= ^iA{Nf + vA{Nf - (3h{Nf - 2aA{N)h{N) + iA{N) - -fb{N) + C 



and the Casimir operator is now realized as 

K = <I>{N + l){(3^-6- 2f3A{N) - AA{Nf)p{N) (2.17) 
+$(A^)(/?2 -6- 2PA{N) - AA{N - l)^)p(A^ - 1) - 2aA{N)\N) 
+ i/3'' -6- 2pA{N))biNf + 2{a/3 - -f)A{N)b{N) + {P^ - 2e)h{N) + ^A{Nf 

+ ^(z/+/x/?)A(iV)3+(-i/./32+^+^+a2+e)A(iV)2+(-i/i/55+^+a7+2C)A(iV) 

Finally the structure function is 

MN) = ) 

p{N-l){AA{N-l)-l){f) + {AA{N) + l){g' 

[{AA{N) + + 2aA{Nfb{N) ~ {[3^ - 6 - 2(3A{N))b{Nf 
-2{a(3- j)AiN)biN) - {(5^ - 2e)b{N) - ^A{NY - ^ (// + pP)A{Nf (2.18) 
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-]^W''-5-2(3A{N)-/\A{Nf){^iA{Nf+pA{Nf-l3h{Nf-2aA{N)h{N)+^A{N) 
with 

f = [5'^ -5-2pA{N)- AA{Nf, g = -5-2(3A{N)- AA{N -If . 

(2.19) 

Thus the structure function depends only on the function p. This function 
can be arbitrarily chosen and does not influence the spectrum. We choose p 
to obtain a structure function that is a polynomial in N, namely we put 

^^^^ " 3 * 2i2/38(Ar + u){l + N + u){l + 2{N + u)f ' ^^'^^^ 

From our expressions for A(N) , b(N) and p(A^), the third relation of the 
cubic associative algebra and the expression of the Casimir operator we find 
the structure function $(A^). For the Case 1 the structure function is a 
polynomial of order 10 in N. The coefficients of this polynomial are 
functions of a, /?, p, 7, 6, e, ly, ^ and C- We give the formula in the 
Appendix. 

Case 2 : For /3 = and 5 7^ we get the solution 

A{N) = V6{N + u),b{N) = -a{N + uy-^{N + u)--^ . (2.21) 



We choose a trivial expression p{N) = 1. The explicit expression of the 
structure function for this case is 

—ae ^ 7^ 7e 07 C ^VS^,^r 
^2(5 A AS 2Sl aVS 2VS 12 



4 4V6 45 25 4 4 
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We use a parafermionic realization in which the parafermionic number 
operator N and the Casimir operator K are diagonaL The basis of this 
representation is the Fock basis for the parafermionic oscillator. The vector 
\k,n >,n = 0,1,2... satisfies the following relations : 

N\k,n >= n\k,n >, K\k,n >= k\k,n > . (2.23) 
The vectors \k,n > are also eigenvectors of the generator A. 

A\k,n >= A{k,n)\k,n >, 
A{k,n) = ^{{n + uy-\-p, f3^0, (2.24) 
A{k,n) = V5{n + u), f] = 0, 5^0 . 

We have the following constraints for the structure function 

<^>{0,u,k) = 0, <l>{p+l,u,k) = . (2.25) 

With these two relations we can find the energy spectrum. Many solutions 
for the system exist. Unitary representations of the deformed parafermionic 
oscillator obey the constraint $(x) > for x=l,2,...,p . 

There are other conditions that should be imposed. The representations 
should be constrained by the differential character of the Hamiltonian and 
the integrals. For example the mean energy should be greater than the 
minimum of the potential 

< H >> minV (2.26) 

This restriction and possibly other ones coming from the differential 
operator character of the integrals should be taken into consideration to 
exclude spurious states. 
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3 Irreducible rational function potentials 



In the case of the three reducible superintegrable potentials the cubic 
algebra is a direct consequence of a simpler algebraic structure. The first 
potential V = '^{x^ + y^) is the well known isotropic harmonic oscillator. 
We can construct the quadratic or the cubic algebra from the Lie algebra as 
in the classical case [10]. The eigenfunctions of the harmonic oscillator are 
well known and are given in terms of the Hermite polynomials. The two 
other reducible potentials V = \{x^ + 2/^) + + p and 

V = ^(4a;^ + y"^) + + cx are two of the four types of potentials found a 
long time ago [6]. There is no Lie algebra in these cases but a quadratic 
algebra [17] and we can obtain the cubic algebra directly from this algebra. 
We obtain from the cubic algebra the same unitary representations that 
were obtained from the quadratic algebra [17]. 

In this section we will apply to the irreducible quantum potentials resuts of 
the Section 2 and give all unitary representations and the corresponding 
energy spectra. Notice that in all case we have f3 = so only Case 2 of 
Section 2 occurs. 

Potential 1 : V = + + r^) 

This potential has the following two integrals 



i2 T-,2 



A = P^-P^ + 2h\^^-^ + -^—^ + ^^^) (3.1) 



+ (3.2) 

^2^^ 4a4 a;2_a2 + (x2-a2)2^'^^• 

The integrals A,B and H give rise to the following cubic algebra and 
Casimir operator 

[A,B] = C, [A,C] = —B (3.3) 

[B,C] = -2h^A'-QteA^H+Sn'H^+Q^A^+S%HA-sK:H^+2-A-2-H- 

a'' a* a* 
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fc4 fc6 fcS fclO 

K = -16h''H'' + 32^H' + 16^H^-40^H-3^ . (3.4) 



The structure function is given by the expression 

^, , ,—ci?E 1,, , .a^E 1-.., .—a^E 3,, , ,—a^E 5,, 

(3.5) 

There are three unitary representations. The first unitary solution is for 
a = ittQ, ao G M. From the condition $(0, u,k) = we find u = — h |. 
The second constraint ^{p + l,u, k) = implies 



2a2 



E = ^ 2 , '^'(a:) = (— )a;(p+ 1 -x)(p + 3-x)(p + 4-x), (3.6) 

where p E N. We have $(p + 1) = which means that the unitary 
representations have dimension p+1. This is also the degeneracy of the 
energy levels. 

The second unitary solution for a = mo, ao G M. We have u = -p n and 



E = -^, $(x) = (-)x(p + l-x)(3-x)(2-x) , (3.7) 

valid only for p=0,l. 
We have 

E > minV = V{0, 0) = — — , (3.8) 

so this is a physically meaningful solution. A third unitary solution exists 
this time for a G M. We have u = + | and 

E= ^^^y , ^x) = {^)x{p+l-x){x + l){x + 3) . (3.9) 



Potential 2 : V = ^{9x^ + y^) 

This potential has the two integrals 

A = P^- + uj\9x^ - y^) (3.10) 
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1 / i2 Q, ,2 

The cubic algebra and Casimir operator of this system are 

[A, B] = C, [A, C] = lAAoj^h^B (3.11) 
[B, C] = -2h^A^ + %h^HA^ - 8h^H^ - ^Quo^ffA + 12uo^h^H 

K = -IQh^H^ + Uuj^h^H^ + 720u;^n^ . (3.12) 

The structure function is 

(3.13) 

4W = (-36.^ft*)(x + „-(^ + i))(x + «-(J^ + l))(. + «-(^ + i)) 

, E 5,, 

We use the two constraints given by the Eq.(2.25). We obtain u = + ^ 
and three unitary representations with the corresponding energy spectra 

2 2 1 

E = 3ujh{p+-), ^x) = {3Quj^h'^)x{p+l-x){p+--x){p+--x) (3.14) 

o o o 

2 4 

E = 3iuh{p + l), (^{x) = {3QLU^h^)x{p+l-x){p+--x){p+--x) (3.15) 

3 3 

4 5 4 

£; = 3c<j%+-), <l>(x) = (36cj2/i4)a;(p+i_a;)(p+__a;)(p+__a;). (3.16) 

3 3 3 

These results coincide with those obtained by solving the Schrodinger 
equation using separation of variable. The eigenfunctions are well known 
and given by 
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where are Hermite polynomials. The corresponding energy spectrum is 

E = ujh{'iki + k2 + 2) (3.19) 
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Potential 3 : F = \{Qx' + y') + ^ 
The two integrals of this potential are 

y 



A = Pl- + io\9x^ - 2/2) ^ (3.20) 



The cubic algebra and the Casimir operator are 

[A, B] = C [A, C] = lUu^ffB (3.21) 
[B, C] = -2n^A^ + Qh^HA^ - Sh^H^ - Su'^h^A + 72u^h^H 

K = -IGh^H^ + 256uj^h'^H^ - lOOSuj^rf . (3.22) 

The structure function is 

(3.23) 

/ / E 7,, 

Using Eq.(2.25) we obtain u = ^ + | and two unitary representations 
$(x) = (36wV)x(p+^-x)(p+l-x)(p+^-x), E = 3uh{p+^), (3.24) 

o o o 

$(x) = (36a;2;i4)x(p+^-a;)(p+|-x)(p+l-x), = 3a;%+l) .(3.25) 

These results are corroborated by those obtained when we use separation of 
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variable and solve the Schrodinger equation. The eigenfunctions are well 
know are given by 

where is a Laguerre polynomial. The corresponding energy spectrum is 

E = ujh{3ki + 2k2 + 4:) (3.28) 



Potential 4 : V = h^-^ + + 

The two integrals are given by the formulas 

,9a;2-y2 i i 

— H \ 

la'^ {y — ay {y + a) 



A = Pl- Pi + 2h\^-^ + J-—- + J-—-) (3.29) 



+ 2 ^^^(y2_a2)2 a^^' 

The cubic algebra and the Casimir operator are 

3fi/?^ 

[A,B] = C, [A,C] = ^B (3.31) 
[B, C] = -2h^A^ - m^A^H - 8h^H^ + 10^ A + 18^H - 24^ 

fc6 fc8 felO 

K = -16n'//^ + 112^/^2 + 9611.// -171^ . (3.32) 

a" ^ ^ 

The structure function is 

$(x) = -^(x+.-(^--))(x+.-(^+-))(x+.-(^+-))(x+.-(^+ 

(3.33) 
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For the case a = iao, oq G M we get the three following unitary 
representations 

$(x) = —x{v + 1 - x){x --){x- -), E = — ^ (3.34) 

QjQ o o ^0 

$(x) = ^x(j> + l-x)(x + ^)(x-i), i? = f^i^±^ (3.35) 

QjQ o o 

<l>(x) = ^x(p + l-x)(x + -)(a;+-), = . 3.36 

For the case a G M we get the three unitary representations 

"^(^) = —r^iP + 1 - 2:)(p + ^ - x)(p + ^ - x), E 

3 3 

'^{x) = —x{p+l-x){p+--x){p---x), E 

3 3 

$(x) = — x(p+l-x)(p+--x)(p---x), £;= .(3.39) 

The Potential bV= h\^[{x^ + y^) + ^ + + j;^] and Potential 6 

V = h^iiiix^ + y^) + + + + (;^] are particular. Their 

integrals of motion A,B and C do not close in a finite cubic algebra. Closure 
at a higher order remains to be investigated. In these cases, we have the 
separation of variables and the unidimensional parts are related to the 
Potential 1 and Potential 4 and their spectra. We will see also in the next 
section that we can obtain information using supersymmetric quantum 
mechanics. 



2a2 



(3.37) 
(3.38) 



4 Supersymmetric quantum mechanics 

In this section we will investigate a relation between SUSYQM [24] and 
superintegrable systems with a third order integral of motion . Let us recall 
some aspects of supersymmetric quantum mechanics. We define two first 
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order operators 

We consider the following two Hamiltonians which are called 
"superpartners" 

(4-2) 

There are two cases. The first is Ail)^^^ ^ 0, E^^ ^ 0, ^Vo^^ 7^ and 
eJ^) ^ 0. We have 

Ef=E«>0, ^i^) = ^M^), V^« = ^AVi^) . (4.3) 

and the two Hamiltonians are isospectral. This case corresponds to broken 
supersymmetry. 

For the second case the supersymmetry is unbroken and we have Aifj^^^ = 0, 
E^i^ = 0, A^i/jf^ ^ and Ej^'^ ^ 0. Without lost of generality we take Hi as 
the potential having the zero energy ground state. We have 



(4.4) 



b2i Je!^' 



^H^' ') Q-(° f) . (4.5) 



We can define the matrices 

Hi 0\ _ 

H2) ^ \A Oj ^ \0 

We get the relations 

[H,Q] = [H,Q^] = 0, {Q,Q} = {Q\Q^} = 0, {Q,Q^} = H . (4.6) 

We have a sl(l|l) superalgebra and Hi and H2 are superpartners. By 
construction, all our potentials can be viewed as the sum of two one 
dimensional potentials H = H^ + Hy. The unidimensional parts of the three 
reducible potentials and the irreducible Potential 2 and Potential 3 are 
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known in SUSYQM. These potentials have the shape invariance property 
[29,31]. We will show that Potentials 1,4,5 and 6 can be also discussed from 
the point of view of supersymmetry. 



4.1 Potential 1 

The Hamiltonian is 

H = H., + H, = ^ + ^ + h\^—-^+ , \, + , \j . (4.7) 



Let us define the two operators 



^ ( h ^ ^ h( ~^ ~^ )) (4 ) 

y/2 dx 2a? x-a x + a 

1 ,^ d h ^/ — 1 —Inn / n 

b = ^ih—-—x-hi + ^— • 4.9 

dx 2a? x-a x + a 

For a = iao, ^ 11^ we have 

' ~ 2 ^ 8a^ + (x - 2ao)2 + (x + zao)^ ^ Aal ^ ^ ^ 

These two unidimensional Hamiltonians are almost isospectral. Hi has a 
zero energy ground state. The supersymmetry is unbroken. This potential 
was discussed in Ref. 33. Non singular superpartners of the harmonic 
oscillator were discussed in Ref. 34 and 35. Coherent states of superpartners 
of the harmonic oscillator have also been studied [36]. Wee see that 
Hi = Hx + ^ IS the Hamiltonian that we are interested in and its 
superpartner H2 corresponds to a harmonic oscillator. 

We apply results for the unbroken supersymmetry. The zero energy ground 
state satisfies 60o = and is 



3 2 1 

b,{x) = al{-)-^-^—^ . (4.12) 



vr 05 + 
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The other eigenfunctions of Hi are obtained by the equation 

(j)^^li = —^=b^(f)n \ In this case ijjn^ are only the eigenfunctions of the 

harmonic oscillator 

{H2) that are written in terms of Hermite polynomials. We get directly for H 



00 , 1 >l 1 -iJ,(3:' + 3™S)„ 2*:, 

-ttfci 7= A/^fci-lJ, 



A=l for ki > 1, A=0 for ki=0. With this expression we get for ki = 

1 -f^a;(3a^ + x2) 

^1 = — e 2 ° 2 • 

V3(27r)5a| «o + 



(4.14) 



We have the following energy spectrum for H 



fc2 

^^^=0, <+i = ^(^i + 3) • (4.15) 



We thus obtain the spectrum of ( the x part of the irreducible Potential 
1): 



^o=-772' ^.^;+i = 7772(^1 + 0) • (4.16) 



If we add Hy to these results we get the energy spectrum and the 
eigenfunctions of the Potential 1. There are two families of solutions. The 
first corresponds to the energies 



{ki + k2 + 2)h^ ip + 2)h' 
^ = 24 = ^ ^ ^ 
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with eigenf unctions 

(4.18) 

and is also obtained from the cubic algebra. The second corresponds to the 
energies 

E = W ' (4.20) 



2al 



with the corresponding eigenfunctions 



3 2 1 6^"^° 

i:{x,y) = 4)o{x)xk2{y), M^)=ao{-)^—^ — ^ (4.21) 

TT + X^ 



and Xk2{y) cis in Eq.(4.19). 

The two states obtained from Eq.(3.8) are given by Eq.(4.20) for A;2=0,l. 
For k^y 3 there are common eigenvalues given by Eq.(4.17) and Eq.(4.20) 
and therefore the degeneracy is p+2. 

Let us consider the case a G M. We have the following Hamiltonians 

Hi = b%=^ + — — + + - — 4.22 

p2 fc2 2 crfc2 

2 8a^ 4a2 ^ ' 

This case is more complicated because of the singularities on the x-axis for 
the Hamiltonian Hi. We have a regular Hamiltonian connected to a 
singular one and we have also for H2 negative energy states. Such situations 
have attracted a lot of attention and many articles were devoted to such 
singular potentials. An important case is the one of Jevicki and Rodrigues 
[37,38]. The corresponding Hamiltonians are 

H_ = ^^+x'-3, F+ = --f- + x2 + 4-l. (4.24) 
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Factorization of Hamiltonians Hi and H2 given by Eq.(4.22) and (4.23) give 
us an algebraic relation that does not take into account the presence of 
singularities or boundary conditions. The wavefunctions given in Eq.(4.3) 
and (4.4) do not necessarily belong to the Hilbert space of square integrable 
functions. The potential in Eq.(4.22) has impenetrable barriers coming from 
the singularities. We can consider the superpartner to be the harmonic 
oscillator with two infinite barriers (at x = ±a) to recover the 
supersymmetry [39]. In the Ref. 39 a superpartner of the harmonic 
oscillator with one singularity was considered but the method can be 
extended to more singularities. The only case that was solved analytically 
and where the energy levels are equidistant is when the singularity was at 
the origin. In our case we were not able to solve analytically and we leave 
for future investigations these numerical calculations that appear 
interesting from a phenomenological point of view. Singular potentials were 
also investigated by A. Das and S.A.Pernice [40] by means of the 
regularization method. M.Znojil [41] has discussed another method that 
consist in the complexification of the potential. In Section 6 we will discuss 
the complexification of the irreducible quantum superintegrable Potential 1. 



4.2 Potential 4 

We apply these results to the next irreducible potential 

V = h\ T^ + i + 7 7:7 4.25 

We can also use SUSYQM because the y part is the same as the x part of 
Potential 1. For the case a = iao, ao G M we find with energy 

E = -^{3h + k2 + 3) (4.26) 



with the corresponding eigenfunctions 



XkAx) = ^^^{7^)--e~^o H,U—^x) (4.27) 
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, , . ao , 1 ,4 1 -^, (?/' + 3?/ag) ^ 2A;2 

(4.28) 



and we get from the singlet state the energies 

2ao 

and eigenfunctions 



E = ^{3ki) (4.29) 



3 2 1 e ""o 

^(x,?/) = 0o(y)Xfci(a^), 0o(?/) = ao(-)^^- — ^ (4-30) 



and Xfci(x) as in Eq.(4.27). 



4.3 Potential 5 

The potential is 

^ = + + / ^9 + / ^ J (4-31) 

For the case a = iao, Oq G M we have 



E=ih±l^ff (4.32) 



with the eigenfunctions given by 



^0 . 1 xi 1 g-f^( (^^ + 3xag) ^ 2h 



V(x) = ^=^(;^)3^==e H{^-^^-L^—^H,,-^XH,,.,) 

(4.33) 

where ^^^(2;) are Laguerre polynomials 
We have also the energies 

2a, 



E = " '^^2 (4-35) 



-^0 
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with the corresponding eigenfunctions 



3 2 1 e^'^i 

ip{x,y) = Xk2{y)Mx)^ M^) = ^oi-)''—^ — ^ (4.36) 
and Xkiiy) as Eq.(4.34). 

4.4 Potential 6 

We consider the potential 

V = + V') + + + + (4.37) 



For the case a = iao, Cq G M we have the energies 



with the eigenfunctions given by 



ao , 1 ,1 1 -^Ax^ + 3xal) 2h 



(4.39) 

«o , 1 ,1 1 -j^, (?/^ + 3j/ag) ^ 2h 

(4.40) 



The singlet state in the x part of the Hamiltonian gives the energies 



nj = 

with eigenfunctions 



3 2 1 e*"o 

^(2^,2/) = Xfe2(?/)'/'o(a^), 0o(2;) = ao'(-)4^- — r (4.42) 

TT an + 
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and XkM as Eq.(4.40). 

We also obtain another kind of solution from the singlet state in the y part. 
The energies are 



h:\k 

with the corresponding eigenfunctions 



3 2 1 e^^o 
TV a^ + y^ 

and (f)ki{x) as Eq.(4.39). 

and a state coming from the singlet state in the two parts with energies 



^0 

with the following expresion for the eigenfunctions 

2 2 

-x^ -y 

3 2 1 ^ 2 1 e*^ 

5 Generating spectrum algebra 

The supersymmetry allows us to find the creation and annihilation 
operators of the x part of the irreducible Potential 1. They are given by 

M = b'^cb, = b^c^b (5.1) 

where c and are annihilation and creation operators of the superpartner 
H2 that is a harmonic oscillator. We have 

and 

^2 dx 2a"^ X — a x + a la^ dx v2 dx la^ x 

(5.3) 
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MT = —(-h- :rX+h( \ )) — (x-2a^—) — (h- :rX+h( 

dx 2a^ x-a x + a 2a? dx\/2 dx 2a? x 

(5.4) 

The zero energy ground state given by the Eq.(4.12) is annihilated be the 
annihilation operator but also by the creation operator. 

M0o(x) = MVo = (5.5) 

The creation and annihilation operator for the y part {Hy) of the Potential 
1 are 

We have the commutators 

[^,^1 = ^(^+^^)^-^(^+^^)-^, [L.L^] = l ■ (5.7) 
We consider the following operators [16] 

E+ = M^L\ = ML, F+ = (M^)^ F_ = G+ = {L'^f, = L 

(5.8) 

We add to these operators the Hamiltonian, the integrals of motion A, B 
and C (Eq.(3.1) (3.2) and (3.3)). We have the following quintic algebra that 
contains 45 relations where the cubic algebra appears as a subalgebra : 



(5.9) 



a^' 



[H, F±] = ±(^)F±, [H, G±] = ±(^)G'±, 



[A,B] = C, [A,C] = —B, [A,E^] = 0, [A, F±] = ±(— )F±, 

a^ 

[A, G±] = t{^)G±, [B, C] = -2h^A^ - Gh^A^H 
a^ 

fc4 fc4 fc4 fc6 fc6 fcS 

a^ a^ a^ a* a'^ a° 



[5, EA = -2^hF. + ^(i/ + \AfG_ - ^-^{H + \a)G_ - 
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[B, F.] = m{H + ^AfE. -^{H+ ^A)E. + ^E., 

[B, = Mh{H + ^-AfE^ - '^{H + \a)E^ - 
[B, G+] = -AihE+, [B, G_] = -AihE_, 

Q„' fc3 

[C, G±] = T^E±, F±] = 0, [i?±, G±] = 0, 

_ 2fc2 o 2fc2 fc4 

^ ' +^ 16 4 8 

Ai^ i^^ H -iJ H -, 

i^-] = ^^C{H + l/l)^ + ^-^B{H + + + \a) 



8a2 ' 2' Ua? 32a4 



8a2 ' 2' 64a2 32a4 
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This polynomial algebra is the spectrum-generating algebra. 

6 Complexification of superintegrable 
potentials 

In quantum mechanics textbooks the Hermiticity of the Hamiltonian is 
often presented as a condition for the energy spectrum to be real. There 
exist other requirements that can be chosen without loosing essential 
features of quantum mechanics. One requirement that appears more 
physical is the space-time reflection symmetry i.e. the Hamiltonian is 
invariant under the PT transformation [42], i.e. the simultaneous reflections 
P : x^-x, p^-p and r : x^x p^-p, i^-i. For potentials invariant under 
such transformations the energy spectrum can also consist of 
complex-conjugate pairs of eigenvalues. The PT-symmetry is thus said to 
be broken. The notion of Pseudo-Hermiticity was introduced by 
A.Mostafazadeh [43]. He shows also that every Hamiltonian with a real 
spectrum is pseudo-Hermitian and that all PT-symmetric Hamiltonians 
studied belong to the class of pseudo-Hermitian Hamiltonian. The 
replacement of the condition that the Hamiltonian is Hermitian by a 
weaker condition allows us to study many new kinds of Hamiltonians that 
would have been excluded and from a phenomenological point of view may 
describe physic phenomena. The case H = + x^iixY was studied in 
detail by C. Bender in 1998 [42]. 

Complexification has been proposed as a natural way to regularize singular 
potentials [41]. It consists in a transformation of the type x — s> x - ie applied 
to a potential. The harmonic oscillator and the Smorodinsky-Winternitz 
potential are PT-symmetric Hamiltonian after a complexification [41]. 
We will consider the complexification of the Hamiltonian 



y ~ 




y {x-ief + {y-ief 
' ^ 8a4 




1 1 
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The complex harmonic oscillator Hamiltonian Hy is known to be 
PT-symmetric. It's energy spectrum is real, namely : 



The eigenfunctions are 



(y-i€)^ ill — if] 

= N^e-^^Hm r ^ ) (6-3) 
v2a 

(here and below Nm is a normalization constant). 

To get the energy spectrum and the eigenfunctions of we complexify the 
operators given by the Eq.(4.8) and (4.9). We get two PT-symmetric 
Hamiltonians. This transformation allows to regularize when a G M. The 
(real) energy levels and eigenfunctions of the Hamiltonian H2 are known. 

The Darboux transformation is still valid for non-Hermitian Hamiltonians 
but supersymmetry is replaced by pseudo-supersymmetry [44]. We have 
b ijjgr = that correspond to the zero energy state of H2 

ifjgr = Ngre~^^^{a^ - {x -ief) (6.6) 

We can obtain the eigenfunction of Hi by applying b' on other state of H2 
given in terms of Hermite polynomials. We get 

^(£q£)!^ 2{x-ie) {x - ie) 

cP. = ^ne-^{j^^--^K.,s{—^) (6.7) 

2(n + 3) Ax-ie) 
Let us give the explicit expression for the ground state and the first excited 
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state 



^o = iVoe^^V^-^ ^ 6.8 



, (3a^ + 2a^i(3: - ie) + {x - ie)^){x - ze) 
(pi = ^ J r^^T 6.9 

The probabilistic interpretation of the wave function of non-Hermitian 
quantum systems [45] is given by a pseudo-norm that is not positive 
definite. 

/oo 
dxip*{-x)^p{x) = a, a = ±1 (6.10) 
-oo 

The corresponding energy spectrum is given by 

E. = 4^ (C.ll) 

We obtain for the complexified superintegrable potential the energy 
spectrum 

with eigenfunction given by Eq.(6.3) and (6.7). 



7 Conclusion 

The main result of this article is that we have constructed a Fock type 
representation for the most general cubic algebra generated by a second 
order and a third order order integral of motion by the means of 
parafermionic algebras. We present in detail the cubic algebra for all 
irreducible quantum superintegrable potentials, the unitary representations 
and the corresponding energy spectra. All cases with finite cubic algebras 
belong to Case 2 of Section 2. Thus they correspond to /5 = in Eq.(2.5) 
and the structure function is given by Eq. (2.22). In two cases of irreducible 
potentials the integrals of motion do not close in a finite dimensional cubic 
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algebra. It could be interesting to see what kind algebraic structure is 
involved in these cases. Comparing with a earlier article [10] we can see 
from this article how the cubic Poisson algebra is deformed into a cubic 
algebra in quantum mechanics. 

The method that we use to find energy spectra with the cubic algebra is 
independant of the choice of coordinate systems. We could apply these 
results in the future to systems with a third order integral that are 
separable in polar, elliptic or parabolic coordinates. The method is also 
independant of the metric and could be applied to superintegrable systems 
in other spaces. The methods developed in this article could be applied to 
other physical systems. One such system is a Schrodinger equation with a 
position dependent mass [32], others arise in the context of supersymmetric 
quantum mechanics. 

The Potential 3 is also a special case of the following potential [46,47] 

V=^{ex' + mV) + ^^+^ (7.1) 

In general this system has integrals of motion of order greater than 3 and 
the more complicated polynomial algebra should be studied. 
All the potentials considered in this article can also be viewed as the sum of 
two one-dimensional potentials H = + Hy. We have investigated each of 
these unidimensional potentials in terms of supersymmetric quantum 
mechanics. The superintegrability of these two-dimensional potentials seems 
to be related to the supersymmetry property. Using the supersymmetry we 
have obtained the energy spectra and the eigenfunctions. We have 
compared the results with those obtained using the cubic algebras. One 
particular feature is the appearance of singlet states. For the Potential 1 
there is an additional degeneracy that is not obtained by the algebraic 
method using the cubic algebra. 

It was shown that many well known potentials such Dirac delta and 
Poschl-Teller display a hidden SUSY where the reflection (parity) operator 
play the role of the grading operator [48]. Potentials with elliptic functions 
can also be discussed from this point of view [49]. Potentials with elliptic 
functions appear in Ref. 9. These cases are not truly superintegrable since 
there exists a syzygy between the Hamiltonian, second order integral and 
the third order integral of motion but it has been shown that the third 
order integral can be used to obtain the eigenfunctions and the spectrum 



29 



[50]. We leave quantum potentials involving Painleve transcendents for a 
future article. 

Superintegrable potentials and their integrals of motion can be complexified 
and investigated from the point of view of PT-symmetric quantum 
mechanics. The complexification appears also as a natural way to regularize 
the singular potentials. 

It's would be interesting to investigate the relation of pseudo-Hermitian 
Hamiltonians and supersymmetry with superintegrable systems. 
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8 Appendix 

Structure function for the case /9 7^ 

(Al) ^{N) = 384/i/3^°A^i°-1920///9^°A^^+(-15365/i/3^+1024i//3^+3040/i/5^° 

-2304/3®a2)A^*+(6144(5^/3^-4096z//3^-640yU/3^°+6144/3^a2)iV7^^2304(5V/3®-3072;/5/?^ 

3072e/3^-76805;u/3^+5120z//3^-2512///3i°-30725/3^a2^2304/3V+3072/5^a7)iV6 

+ (-69125V/3^+9216z/5/5^-9216C/9^+1536(5/i/3^-1024///3'^+1712///9i°+92165/3^a2 

-7680/5^a^-9216/3^a7)iV^+(-15365V/3^+3072z/52/?^-6144e5/9^+63365V/5'' 
-8448z^5/3^a7 + 8448^/3^ + 32646^(3^ - 2176z//9^ + 428///?^° 
+A6085^/3^a^ - 8U86p^a^ + Q72(3^a^ - 92165/3^07 + 8448/3^ 
+3072/3^2 + 6144/5^ae + 12288/3^C)A^^ + (30725^/5^ - QlUuS^P^ 
+ 12288^(^/5^-1 1525 1536i/(5/?^- 1536^/?^- 19205/^/3^+ 1280i^ (3^ -616 ^f3^^ -122886'^ (3W 
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+1536(5/3^a2+2688/5^a2^24576(5/?^a7-1536/3^a7-6144/?V-24576;5^ae-24576/3^C)^^^ 

119 

+ 1728u6(3'^-966fi(3^+64u(3^+—^f3^^-12288f3^K-3072S''^(3^a'^+69126^f3^a^+1728S 

/3''a2-624/3^a2^92165^/3^a7-138245/3^a7-1728/9^a7-61445/3V+1536/?V 
-122885/9^ae+9216/3^ae+12288/3V-12288(5/3^C+18432/?^C)A^^+(-3845V/5^ 
+1024u6^f3^ - 3072^5^/3^ + 2565^/?^ - 512u5^(3^ + 3072^5/3^ + 2085^/9^ 

129 

-960z/5/9^+960e/5^+2885^/5^-192z//9^+— ///9^°12288/9^/r+30725=^/32a2_9go5/5^a2_ 

288/3^a2 - 92166'^ P^aj + 960/3^a7 + 61445/?V + 1536/3Vl22885/3^ae 
+6144/3'^ae-12288/3V+122885/9^C-6144/3^C)^^+(965V/3^-256z/(5^/5^+768<e5^/?^ 

+325^/3^ - 64//52/35 - 384^(5/3'^ + 206^fif3^ + lUud(3'^ - 144^/5^ - 546 j^P^ + 36i^(3^ 
117 

H(3^^ - 3072/5^ + 768<5^a2 _ ygs^^^^^^ _ 430^2^4^2 ^ 144^^6^ + 87P^a^ 

8 

-30726^ (3a-f + 23046^P^a-f + 9Q0S(3^a-f - UAp'^a-/ + 30726^ P^-/^ 
-15366 P^-f^ - 576/?%^ + QlU6'^p^ae - 30726 pW - 2688/?^a 
e - 122886p^^e + 3072p^-fe + 12288/3^^ - 3072(5/?^C + 768P\ 
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